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For very general scalar-field theories in which the equations of motion are at second-order, we 
evaluate the three-point correlation function of primordial scalar perturbations generated during 
inflation. We show that the shape of non-Gaussianities is well approximated by the equilateral 
type. The equilateral non-linear parameter f^ 1 is derived on the quasi de Sitter background 
where the slow-variation parameters are much smaller than unity. We apply our formula for f^ 1 
to a number of single-field models of inflation-such as k-inflation, k-inflation with Galileon terms, 
potential-driven Galileon inflation, nonminimal coupling models (including field-derivative coupling 
models), and Gauss-Bonnet gravity. 

I. INTRODUCTION 

The inflationary paradigm ll| can successfully account for the observed temperature anisotropics in Cosmic Mi- 
crowave Background (CMB) [2| as well as the galaxy power spectrum [3|. Although we have not identified the origin 
of inflation yet, the upcoming observations such as PLANCK [4] are expected to provide more high-precision data for 
_ ■ the primordial scalar/tensor power spectra. This will allow us to distinguish between a host of theoretical models of 
!jj inflation [fj| observationally. 

^O ■ In addition to the spectral index of scalar density perturbations n-jz and the tensor-to-scalar ratio r, the information 

of primordial scalar non-Gaussianities is useful to break the degeneracy between different models (a-Q ■ I n standard 

slow-roll inflation in which cosmic acceleration is driven by the potential energy of a scalar field (f>, the non-linear 

parameter f^"' characterizing the equilateral shape of non-Gaussianities is small (|/nl U1 I ^ 1) [lOL [FL| (see also 

<~ "i, Refs. |L2| - tl4| ). Meanwhile, in the presence of a non- linear field kinetic ener gy i n X = — 9 M 09 M 0/2 such as k-inflation 

Si' [H|, Dirac-Born-Infeld (DBI) models [MEl, (dilatonic) ghost condensate flUllij . and Galileon inflation (HEH], ii; 

is possible to realize large non-Gaussianities detectable in future observations (I/niI" I ^ 10) [22b25j . 
£jy The viable theoretical models of inflation are usually constructed to keep the equations of motion at second-order 

in order to avoid the Ostrogradski's instability [26]. In 1974 Horndeski [27| derived the most general single-field 
Lagrangian giving rise to second-order field equations. Recently there has been renewed interest in second-order 
\^J gravitational theories in connection with the Dvali-Gabadadze-Porrati braneworld [28[ and Galileon gravity |29l . 1301 ] 
(see also Refs. |3ll - l4G| ). Deffayet et al. [41] showed that the most general action in those theories is given by 
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-^R+ P(</>, X) - G 3 (0, X) □</> + £ 4 + C 5 



(1) 



Here g is the determinant of the metric g^ v , M p \ is the reduced Planck mass, R is a Ricci scalar, and 

U = G 4 (ct>,X)R+G A ,x[(n ( t>) 2 -(V IJ ,V v <f>)(Vi 1 V v <£)], (2) 

£ 5 = G 5 (^X)G^V^V^)-jj;G 5 ,x[(n<£) 3 -3(n<£)^^ (3) 



where Gi (i — 3, 4, 5) are functions in terms of </> and X = —d^(f)d ll (j)/2 with the partial derivatives G^x = dGi/dX, 
and G^ = R^ u — Ojj^R/2 is the Einstein tensor {R^ v is the Ricci tensor). The action (JXJ) is equivalent to that derived 
by Hordenski |42ll43l|. In Ref. J4j] the spectra of scalar and tensor perturbations were derived for the above models. 

In this paper we compute the three-point correlation function of scalar metric perturbations for the inflationary 
models described by the action ([T]). Note that the calculation of scalar non-Gaussianities up to the term Gz(4>, X)E\(f) 
was carried out in Refs. |44TJ46l ] (see also Ref. [47J for non-minimal coupling models with F{4>)R). Taking into account 
the Lagrangians £4 and £5, we can cover a wide variety of gravitational theories such as scalar-tensor theories [48j 
(G 4 = F{4>)), field derivative couplings with gravity [MlHSl (G 5 = G(<f>)), Galileon gravity [HIIj (G 3 oc X, G 4 oc X 2 , 
G5 oc X 2 ), higher-curvature gravity [51] (including Gauss-Bonnet and f(R) theories), a' corrections appearing in low- 
energy effective string theory [52j |. 

Even in the presence of the new contributions coming from the terms G4 and G5 we shall show that the shape of 
non-Gaussianities can be well approximated by the equilateral type. The equilateral non-linear parameter /ni!" ^ s 
derived on the quasi de Sitter background. We also apply our formula to concrete inflationary models in order to 
present the cases in which large non-Gaussianities can be realized. 



This paper is organized as follows. In Sec. [IT] we derive the background equations of motion on the flat Friedmann- 
Lemaitre- Robertson- Walker (FLRW) background and introduce a number of "slow- variation" parameters convenient 
for the calculation of primordial non-Gaussianities. In Sec. Mil the action (JJ is expanded at second-order in pertur- 
bations in order to obtain the power spectrum of curvature perturbations generated during inflation. We also obtain 
the spectrum of gravitational waves and the tensor-to-scalar ratio. In Sec. HVI we expand the action (TTJ) at third-order 
in perturbations and compute the three-point correlation function of curvature perturbations. In Sec. [V] we study the 
shape of non-Gaussianities for the new terms appearing in the three-point correlation function. In Sec. IVII we derive 
an explicit formula for the equilateral non-linear parameter /^" under the slow- variation approximation. In Sec. IVHI 
our formula for J^l" i s applied to a number of concrete models of inflation. Sec lVHil is devoted to conclusions. 

II. BACKGROUND EQUATIONS 

We first derive the background equations for the theories given by the action (fl) on the flat FLRW background 
with the scale factor a(t), where t is cosmic time. Taking variation of the action at first order with respect to the 
metric elements <?00i 9u, an d the field 0, where goo corresponds to the Lapse function N , we find 

E l = 3M 2 l H 2 F + P + 6HG i ^+(G 3 ^-12H 2 G^ x + 9H 2 G 5 , lj) -P x )ft 

+ (6 G^x - 3 G 3 , x - 5 G^ X H 2 ) Hft + 3 (G 5 , x - 2 G i>xx ) H 2 ft - H 3 G 5 , X xft = , (4) 

E 2 = 2[(G 5 ^-2G^ x )ft-HG 5 , x ft + FM 2 l }H + 3M 2 l H 2 F + P + 4HG^<j> 

+ [2 G 4 , + Lff (G 5 , - G^ X )<P + (2 G^x - G 3 ,x - 3 H 2 G b , x )ft + 2H{G 5 , <t>x - 2 G^ X x)ft - H 2 G^ X xft]4> 
+ (2 G 4>0 + 3 ff 2 G 5 , - G 3 , - 6 H 2 G A ,x) ft + 2H (G 5 , 00 - G b , x H 2 - 2 G^ x ) ft - H 2 G 5 ^ x ft = , (5) 

E 3 = [6G 4> + 12(G M - Gi,x)H4> + 3(2 G^ x - G 3 , x - 3H 2 G 5 , x )ft + 6(G 5j0 x - 2 Gi. X x)Hft - ZH 2 G b , X xft]H 
+ {3(G5 : 0xx — 2Gi^xxx)H 4> — H G$^xxx<t> + [3(2G 4:< /,;or — G 3 ,xx) — 7 H G^^xxjH^ + 2G 3 ^ — ?j 
+ [3(5 G^x - 8 G 4 .xx)H 2 + G 3 .^ x - P.xx]ft + 6(3 G^ x - G 5 , X H 2 - G 3 , X )H4> + 6(G 5 ,^ - G^ X )H 2 }4> 
+ 3(G5^^x — H Gs t xx — 2Gi^xx)H 4> + [KjGn^x — 1&G4,,xx)H + 3(2G 4j 00jf — G 3 ^x)]H(f> 
+ [3(G5,00 + 6 G 4 ,0x — 3 G 3 .x)H —9H G$^x — P,<px + G 3 ^](f) 
+ 3[6(G 5 , - G^x)H 2 + 2 G 3i0 - P,x]HJ> - H 3 G 5 ^ X xft + 12 H 2 G^ + Pj = , (6) 

where H = a/a is the Hubble parameter, and 



2G 4 



(7) 



Note that a dot represents a derivative with respect to t, whereas a comma corresponds to a derivative in terms of 
(j) or X (e.g.. G^^x = d 2 G^/dXd4>). These equations, because of the Bianchi identities, are not independent as one 
can directly verify that 

4>E 3 -E x - iH{E 1 -E 2 ) = 0. (8) 

Eliminating the terms P from Eqs. ^ and (|5J), it follows that 

(1 - 4<5 G4 x - 2S G5X + 2S G5t p)e = 5 PX + 3S G3 x ~ 25 G3 <t> + 6 8 GAX ~ S Gil p - 6 S G5l p + 3 S G5X + 12 Sgaxx + 2 S G5X x 

— 10 Sa^x + 2 5c4<t><t> — 8 5g5<px + 2 5c5<t><}> — S^(5g 3 x + 4 Sgax — <>G4if> 
+8 Sgaxx + 3 Sg5X - 4 S G 5<p + 2 $G5XX - 25ga<px ~ 4 o~co<px) , (9) 

where we have defined the slow- variation parameters 

H (j) x P.xX G 3> x<f>X G 3 . ( j > X _ d.xX 

z = -jp> ^'ITft 5px ~1^pf' 5g ' x '1^hf' 6g3 *-m^ipf' fc4x= Mp r ' 

G 4 ,0</> G4. c f,x4'X G^^X G±.xxX G^^X 

^-^HF-' 6GA * X ~ M 2 ,HF ' SG4 ^ = M^IPF' Sg4XX ~ M 2 pl F ' 5 °^ = l^F ' 

r G^,.xH(j)X G5 : xxH<t>X G^^xX G^ t ^<pX 

6G5X = M a p , 6 G5 xx- jjp 5 GB 4,x= M y , 6gsm> = M ^ HF ■ ( 10 ) 



Since e«l during inflation, we require that the slow-variation parameters defined above are much smaller than the 
order of unity. 

Taking the time-derivative of the quantity 5gax , we obtain the first-order quantity 

Sgax 26G4xx8cf> , 5ci(t,x . or £ ,-.-.^ 

VG4X = 77T = — 7 + -t h 2<^ - S F , (11) 

U0G4X 0G4X OGAX 

where 6f = F/(HF). This means that 6g4<j>x is higher than the first order. Likewise one can show that 

{<)G3tj)X, ^G300, <5g40X, ^G4<^, $G5<pX, SG5<p<p} = 0(e ) , (12) 

where 



r _ G^^xX _ G3., 

G30 * _ M 2 H 2 F ' G3<W ~ M 2 H 3 F ■ 



(13) 



'pi- 1 - 1 - 1 J "pi J 

Then, at first order, Eq. © reduces to 

e ~ 5p X + 38 G 3X - 2<5 G 3</> + 6 8 G ax - <5g4</> ~ 6 (5 G5 + 3 S G 5X + 12 fexx + 2 <5 G5 xx • (14) 

It is also convenient to notice the following relation 

S F = 2<5 G40 + tecAxSj , (15) 

that we will use hereafter. 

III. SECOND-ORDER ACTION 

In order to discuss the primordial non-Gaussianities we need to first study linear perturbation theory. We shall 
use the momentum and Hamiltonian constraints to integrate out all the auxiliary fields. For the calculation of scalar 
non-Gaussianities it is convenient to choose the ADM metric [53| in the form 

ds 2 = -[(1 + a) 2 - a{t)- 2 e- 2n (dip) 2 } dt 2 + 23^ dt dx i + a(t) 2 e 2K dx 2 , (16) 

where a, tp, and 1Z are scalar metric perturbations [Il|. We choose the uniform field gauge, Sip — 0, to fix the 
time-component of a gauge-transformation vector f^ . The spatial part of £ M is fixed by gauging away a perturbation 
E that appears as a form Eq in the metric (|16p [54| . 

Expanding the action ([IJ up to second order in the perturbations, we find the following result 



S 2 = / dtd 6 x a 6 
where 



-3wi7t 2 + -^ (2wiR, - w 2 a)d 2 ip - -^- ad 2 K + 3w 2 aTZ + - w 3 a 2 + -^ (dllf 



(17) 



tui = M 2 yl F - AXG^x - 2HX<}>G 5tX + 2AG 5> , (18) 

w 2 = 2M 2 l HF-2X^x-16H(XG i ,x + X 2 G^xx) + 2^(G i ^ + 2XG^x) 

- 2H 2 <p(5XG 5 , x + 2X 2 G 5 ,xx) + 4ff A(3G 5 , + 2XG 5)0 x) , (19) 

w 3 = -9M 2 ,H 2 F + 3{XP X + 2X 2 P XX ) + 18HJ>(2XG 3)X + X 2 G 3 , XX ) - 6A(G 3 , + XG 3 ,*x) 
+18H 2 (7XGi. x + l6X 2 G^xx + 4X 3 G^ X xx) ~ 18F0(G 4 , + hXG^x + 2X 2 G^ X x) 
+ 6H 3 <fi(15XG 5 ,x + UX 2 G^xx + 2X 3 G, 5X xx) - 1&H 2 X(6G 5 ^ + 9XG 5 ^ X + 2X 2 G 5 ^ XX ) , (20) 

w A = M 2 1 F-2XG 5i(j> -2XG 5i x4>. (21) 

Although the coefficients are quite complicated, the structure of the action (fT7)) is similar to that found in Ref. [45j ■ 
It is straightforward to find the momentum and Hamiltonian constraints 

a = LiR, (22) 

-VV= ^a + m- 2 -^±d 2 K, (23) 



where 

U = ^ . (24) 

Substituting Eq. (|22)l into Eq. (JT7J) and making integrations by parts, the second-order action reduces to the following 
form 



S 2 = I <!!<>' 
where 



n 2 - ^ [&Rf 

a z 



(25) 



wi(AwiW3 + 9w?,) 



<>> = J.2 - ' (26) 

— W2W4 + Aw\w 
Wi(4wiW3 + 9w?,) 



3w^ 
3(2w\w2H — W2W4 + 4wiWiW2 — 21V1W2) 



ci = ^^ ^^, _ (27) 



In order to avoid the appearance of ghosts and Laplacian instabilities we require the conditions 

Q>0, c 2 >0. (28) 

For later convenience we introduce the following parameter 

_ Qc 2 s _ 2w 2 w 2 H — w 2 Wi + 4wiibiW2 — 2w 2 w 2 , ■. 

Expansion of e s in terms of the slow- variation parameters leads to 

e s = e + 5qzx + 5g44> + 8S GiXX + S G5X + 2S G5X x + 0(e 2 ) 

= 6 PX + 45 G3X - 2<5 G30 + 6S GiX + 20S G4XX + 45 G5X + 4S G5XX - 65 G54> + 0(e 2 ) , (30) 

where we have used Eqs. ([15]). ([II]). and (fT5)) . 

The two-point correlation function of curvature perturbations 1Z can be derived by employing the standard method 
of quantizing the fields on a de Sitter background [54|. The power spectrum of 1Z, some time after the Hubble radius 
crossing, is given by 

H 2 H 2 

Vn = 8^f = WM^FtsCs ' (31) 

Its spectral index, evaluated at c s k — aH (k is a comoving wave number), can be derived as follows 

_ dlnP n 

dlnk c s k=aH 

= -2e-5 F -T ls -s = -2e s - Vs - s + 2S G3X + 16S G4XX + 2S G5X + 4S G5XX + 0{e 2 ) , (32) 

where 



Vs = lk s > S = ^T (33) 

Here we have assumed that c s is a slowly varying function, such that \s\ <C 1. 

Let us now proceed to the power spectrum of the gravitational waves for the theories under consideration. We 
study the tensor perturbations of the form 

ds 2 = -dt 2 + a 2 (t) (Sij + h T J) dx l dtf , (34) 

where hf T is transverse and traceless. It is known that the hf T can be decomposed into two independent polarization 
modes, namely hf T = h + e[: + h x e^. We choose the normalization for the two matrices such that, in Fourier 



space, eh{Je)eh(— k)* — 2. (where A = +, x), and efAk)e*A~k)* = 0. In this case the second-order action for the 
gravitational waves can be written as 



S = Y^ I dtd 3 xa 3 Qq 



h\-% (dh x f 



(35) 



where 

1IU l 

M 2 F (1 - 4<W - 26 G5X + 25 G ^) , (36) 



4 4 pi 



4 = ^ i = l + 4<5 G 4x + 2(5 G 5 X -4 ( 5 G50 + O(e 2 ). (37) 

Provided that {|#G4x|j |<^G5x| 5 |fcs^|} "C 1, one has Qt — M 2 ]F/& and hence the no-ghost condition Qt > is 
satisfied for F > 0. Since c T is close to 1, there are no Laplacian instabilities for tensor perturbations. 
The spectrum of tensor perturbations is given by 

T 2tt 2 Qt4 ~ ^M 2 Y F ' { ' 

together with the spectral index 

d\nV T 



1%T = 



d\nk 



= -2e -6 F = -2e s + 25 G3X + 16S G4XX + 26 G5X + 4S G5XX + 0{e 2 ) . (39) 

For those times during inflation when both Vt and V-r remain approximately constant, the tensor-to-scalar ratio can 
be estimated as 

r = — ~ 16c s e s . (40) 

Then we have the consistency relation 

r ~ 8c s {-n T + 26 G3X + l65 GiXX + 25 G5X + AS G5XX ) . (41) 

The consistency relation in standard inflation (r = — 8c s 7t.t) is modified because of the presence of the terms Gi 
(t = 3,4,5). 

The spectra of scalar and tensor perturbations coincide with those derived in Ref. [43j . 

IV. THREE-POINT CORRELATION FUNCTION 

We proceed to the calculation of the three-point correlation function of curvature perturbations for the theories 
described by the action ([IJ. In doing so we expand this action up to third order in perturbations [11]. Although the 
calculation is quite involved, the steps for the derivation of the three-point correlation function are similar to those 
given in detail in Ref. [45||. The third-order action is given by 



/ dt d 3 x a 3 {oi a 3 + a 2 (a 2 K + a 3 H + a 4 d 2 K/a 2 + a 5 d 2 ip/a 2 



+ a [a 6 diRdiip/a 2 + a 7 TZTZ + a 8 lid 2 K/a 2 + a 9 (didjipdidjip - d 2 ^d 2 ip)/a A 

+ aioidid^didjfi - d 2 iJjd 2 K)/a A + a n K d 2 ip/a 2 + a 12 K d 2 tp/a 2 + a 13 K d 2 K/a 2 + a u {dllf/a 2 + a 15 K 2 ] 

+ bi il 3 + b 2 K {dUf/a 2 + b 3 U 2 K + Cl UdiTZd^/a 2 + c 2 li 2 d 2 ip/a 2 + c 3 KKd 2 ip/a 2 

+ (dill + d 2 TZ) (didji/jdidjip - d^d^/a 4 + d^TZd^d 2 ^/^} , (42) 



where 

ai = 3M 2 ,H 2 F - {XP X + AX 2 P XX + 4X 3 P xxx /3) - 2H4>(10XG 3jX + IIX 2 G 3 , XX + 2X 3 G 3 , XXX ) 

+2X(G 3 . c / > + 7XG 3 ^ x /3 + 2X G 3 ^ xx /3) — 2H ("&2>XG^ X + 126X G^ xx + 68X G^^ xxx + 8X G^ xxxx ) 
+2H0(3G^ + 21XG^ X + 24X 2 G^ XX + 4X 3 G 4 , 4lXXX ) - H 3 <j)(70XG 5 , x + 98X 2 G 5 , XX + 32X 3 G^ XXX 

+8X A G^ xxxx /3) + 2H 2 X(30G 5 ^ + 75XG 5 ^ X + 36X 2 G 5 ^ XX + 4X 3 G 5 ^ XXX ) , (43) 

a 2 = w 3 , (44) 

a 3 = -3a 5 

= -6M 2 { HF + 60(2XG 3 ,x + X 2 G 3 , XX ) + 12H(7XG^ X + 16X 2 G^ XX + AX 3 G^ XXX ) 
-H{G^ + 5XG^ X + 2X 2 G^ XX ) + 6H 2 <j>(l5XG 5 , x + 13X 2 G^ XX + 2X 3 G^ XXX ) 

-12HX(6G 6t f + 9XG 5 ^ X + 2X 2 G 5 ^ XX ) , (45) 

a 4 = -4(IG 4 ,x + 2X 2 G^ XX ) - 8H<p(XG 5 , x + X 2 G^ xx /2) + 4X(G 5 , + 2XG 5 , 4>X ) , (46) 

a 6 = -a 7 /9 = an = -w 2 , (47) 

a 8 = 2a w = 2bi = -2c 2 = -4di = 4X0G 5 ,x , (48) 

ag = Oi2/4 = -015/6 

= -M^F/2 + 4(IG 4 j + X 2 G^ XX ) + H^(5XG 5 , X + 2X 2 G 5 , XX ) - X(3G 5 , + 2XG 5 ^ X ) , (49) 

13 = 2a 14 = 263/9 = -ci = -c 3 = -4d 2 /3 = d 3 = -2wi , (50) 

b 2 = w±. (51) 

We use Eq. (l2"2"j) to eliminate the field a from the action (|4"2"|) . which gives 

53 = / <ft d 3 x ^{Ai^ 3 + A 2 K 2 d 2r Jl/a 2 + A 3 K 2 d 2 ijj/a 2 + A 4 WJZ 2 + (A 5 TZ + A 6 K) (d l d j ^d l d J ^ - d 2 ^d 2 i;)/a 4 

+ A-rfydidjifrdidjR. - d 2 V<9 2 ft)/a 4 + A 8 TZ(dTZ) 2 /a 2 + A 9 d l TZd l ^ d 2 ^/^} , (52) 

where 

Ai = 61 + L1O15 + 1*103 + L\ai , A 2 — L\ (Lia 4 + a 8 ) , A 3 = c 2 + L1O12 + I-iOs , 



Aa = 63 + Lia 7 + L x a 2 = 3Q , A 5 — Lia 9 + di , A 6 = d 



2 



A 7 = L iaw , A s = b 2 + a 13 L 1 /2 + L 1 {d 13 + Ha 13 )/2, A g = d 3 . (53) 

The next step is to eliminate the field tp by using Eq. (|23[) . Introducing an auxiliary field X satisfying the relation 

i/> = -Lift + — , (54) 

Wl 

it follows that d 2 X = QTZ from Eq. (|23J|. Plugging Eq. (O into the action ([52"]). we obtain 

dtd 3 x(a 3 f 1 +af 2 + f 3 /a) , (.'»:,, 

where 

\ W\ w( / \ w( ) w( 

+ -4 (%ft + A 6 ft) (didjX)[did 3 X) , (56) 

2LiQ ^ g^ - 2 „ 2 _ , . 2LiQ • 2 «/-»r,x2 , I '1 ( 3^/q^a2 



a 2 - A3L1 + a 5 — — - a 7 — ft^a^ft + ^ 6 — —nnd z n + A 8 ft(5ft)^ - a 9 — —n(dn) 
+ A? - 2A5Ll -R.(d i d j n)(d i d i x) - ^hh. n(d l d 1 n){d l d 1 x) - ^± d 2 nd l nd l x , (57) 

W\ W\ W\ 

f 3 = (A 5 L 2 - A7L1) K [(didjtyididjll) - {d 2 TZ) 2 } + A 6 L 2 TZ [(d t d J TZ){d l d J TZ) - {d 2r R) 2 ] + A 9 L 2 {dTZ) 2 d 2 TZ .(58) 



Along the same lines of Ref. [45|, [55| we express the action (|55l) in a more simple form by carrying out many 
integrations by parts. Finally we reach the following action 



S 3 = dt£ 3 , 



(59) 



where 
£3 



d A x< a^dM^TZTZ 2 + aC 2 M 2 {JZ{dTZ) 2 + a 3 C 3 M pl 7^ + a 6 C^l(p{R) (d % X) + a A (C^ j M 2 ^d 2 1l{d X) 2 

aC 6 Tl 2 d 2 TZ + C 7 [d 2 TZ(dTZ) 2 - TZd i d :i (d l TZ)(d j TZ)] /a + a(C 8 /M pl ) [d 2 TZdiTZdiX - TZd l d j {d i TZ)(d :i X) 
5C 2 



Fi 



sn 



and the dimensionless coefficients Ci (i = 1, ■ 



Cl - ITT 



M. 



C 2 = 



C, 



pi 



K 



3Q + qi(Q + 3HQ)-Qq! 

1 d 



A 



dt 



(0L1Q) 



_L ( Al + A 3 Q - qi ( 
M p i V w\ 






— (Ae+Ag)-^-^- ( ^ 

dt \Wi 



Wl 



C B = 



M 



pi 



A fi 



d /A 

dt 



d fA 5 

1 

3HA 5 



3HA 5 

Wl 



A 2 - A 3 L X 

Qc 2 



Ci = q 3 -^- JL (A 7 -2A 5 L 1 ). 

ZWi 



Cs 



M, 



pi 



g 2 2c 2 A 5 Q 



The terms qi, 52, and 53 appear during the step to derive (|59|) from (|55[) . and are given by [4£ 

92 : 



ii 

„2 



93 



4A 6 ii 2 d f A 7 - 2A 5 L 1 
w\ dt \ a 2 w\ 

2 ad f A 5 L 2 — A 7 L\ 



2A 9 L 1 



U'l 



A 6 Lf 



3 dt \ a 

The last term in Eq. ([50|) is the product of the following quantities 



A 9 L 2 



(60) 

(61) 

(62) 

(63) 

(64) 

(65) 
(66) 

(67) 

(68) 

(69) 

(70) 
(71) 



Ti = ^| {(d k K)(d k X) - d-tdidjftdiKXdjX)]} + qi KK + Al A 2A f x 



and 



SC 2 

sn 



dt 



(a^QK) - aQcid 2 K 



{{dnf - d^didjKdiiz^djn)}} , (72) 



(73) 



The coefficient T\ includes only the time and spatial derivatives of 1Z and X. In the paper of Maldacena [TT( the 
term 7Z 2 is also present in the expression of T\ , which gives rise to the contribution of the order of the slow- variation 
parameter in the final expression of the non-linear parameter. After integrations by parts it is possible to move this 
contribution to other terms such as C\ [45J|, which we have done in Eq. (|60|. Hence, as in Refs. [56j, we can neglect 
the contribution coming from the term (72). 



Defining the conformal time as r = f a 1 dt, the vacuum expectation value of curvature perturbations for the 
three-point operator at r = 77 is (ill . |24L157| 



(K(ki)n(k 2 )K{k 3 )) 



T f 



dro(0|[^(r/,fci)^(r/,fc2)^(T/,fc3),^i„t(r)]|0), 



(74) 



where the interacting Hamiltonian "Hint is given by "Hint = —£3- Note that r, corresponds to the initial time at which 
the perturbations are deep inside the Hubble radius. One can take r, — > —00 and Tf —> because r ~ —l/(aH) during 
inflation. When the integral with respect to r is carried out, we treat the terms Ci (i = 1, • • • , 8) as constants because 
their variations are small relative to the scale factor a. Each contribution of the three-point correlation function 
coming from the terms C* (i = 1, ■ • • , 8) is |45| 



.CiM^H* 



16Q 3 cf (fcifefe) 3 V K 
C 2 M^H i l 



(ft(fci)ft(fc 2 )ft(fc 3 )} (1) = (2 7 r) 3 5^(fc 1 + fc 2 + fee 
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fc?(fc 2 ■ fe 3 ) 
A 



A 
fci 
A 



+ sym 
sym 
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kik 2 + k 2 k 3 + k 3 ki 3kik 2 k 3 
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8Q 3 c s 10 (hk 2 k 3 ) 3 A 
x [k 2 (k 2 ■ fe 3 ) - (fci • fe 2 )(fei • fe 3 ) + sym] , 

1 1 



A 2 



A 3 



(75) 

(76) 
(77) 
(78) 
(79) 
(80) 

(81) 



C 8 H 5 



32Q 2 Af pl cf (fcifc 2 fc 3 ) 3 A 



2h+k 2 + k 3 2ki(k 2 + k 3 ) 



A 



A 2 



where A = k\ + fc 2 + fc 3 . The symbol "sym" corresponds to the symmetric terms with respect to k\, k 2 , and k 3 . 
We express the three-point correlation function in the form 

(ft(fei)ft(fc 2 )ft(fc 3 )) = (27r) 3 S^(k 1 +k 2 + k 3 ){VufB n {ku k 2 , k 3 ) . 

where V-r is given by Eq. pip , and 

(27T) 4 



[k 2 (k 2 ■ k 3 ) - (fci • fe 2 )(fci • fc 3 )] + sym I , (82) 



(83) 



B K (ki,k 2 ,k 3 ) 



nL*f 



-4tj 



(84) 



Collecting all the terms in Eqs. (|75| -(|82 p we have 
A n = 



"*-p\ J 1 / 2 v~^ ,9,9 1 V> , 9 , 3 \ ,, 1 /I V"^ , 1 . 2 V^ ,9,9 I 
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-Co 



pi \ i iyy ijkj j 
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i*3 
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»^J 



(85) 



We also define the non- linear parameter /nl, as |2l.[58l| 

/nl = 



10 A 



7? 



3 Etr *f 
If the three-point correlation function is described by the equilateral configuration (k\ = k 2 = fe), it follows that 

40 -W-, I IV 1/7 1 C) I / <) \~ 

/nl 



(86) 



40M p 2 



9 Q 
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12 96c : 



17 r +l-^-r .ilr -if-^-i r > _i_i r 

2 f 72 M Dl 3 24 M 2 U 24 I M 2 j ° 5 + 36c 2 ^M pl/ ' 6 



17 H 



13 f H V 

V/ iyZC s jw pl JW pl 



P 1 / 



(87) 



This is the same form as that derived in Ref. [45], but the coefficients Ci ({=!,•■■ ,8) are different. 



THE SHAPES OF NON-GAUSSIANITIES 



The shapes of non-Gaussianities have been discussed by a number of authors (see e.g., Refs. [2J, |59j). In Ref. [24j 
Chen et al. studied the case of k-inflation, where only the terms proportional to Ci (with 1 < i < 5) are present. Here 
we wish to discuss the shapes coming from the remaining terms, that is, we focus our attention only on the following 
terms 



g(6) _ 



1 



m 2 (hkzhY 



E U kf Qcl K " 



c& , 



B. 



(7) 



H 2 1 



Ef=i kf 2Q4 K\ K 2 ^ 



— 2_^ kikj 
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3kxk 2 k 3 
K 3 
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B. 



(8) 

■R 
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E Li k? 8c 2 Mp! K 2 



| fc!fc 2 fc 3 E fc2 - ^ ^^ 2 - 6 E fe * 2fe i - E fc " 



k K Y,* 



(89) 



(90) 



Although the total bispectrum is the sum of all the previous contributions (in addition to the standard ones of 

k-inflation), we study the momentum dependence of each single contribution Bj^ (with 6 < i < 8). In fact, we will 

check that the non-Gaussianities associated to each single term Bj^ can be well measured by means of the equilateral 
estimator defined as 



B^ il = (2tt) 4 ( ^fZ il 



1 



1 



1 



1 



jL3l,3 £.31,3 £-3)1,3 £,2jl2l2 j. 7,21,3 



(5 perms.) 



(91) 



10 




k,/k 




Mk 



Figure 1: The shape functions B$(l, k 2 /ki, k 3 /ki)(k 2 /ki) 2 {k 3 /kif (left) and B%\l, ka/ki, k 3 /ki)(k 2 /ki) 2 (k 3 /kif (right). 
The two functions are plotted in the domain 1 — k 2 /ki < k 3 /k\ < k 2 /ki. The lower boundary is given by the triangular 
inequality, whereas the upper boundary is chosen in order not to repeat twice the same physical configuration. Finally, the 
functions Bj^ are multiplied by the measure (k 2 / 'fci) 2 '(k 3 / 'k\) 2 following Ref. [59l |. The plots are normalized to have a unit value 
at the point k 2 jk\ = k 3 jk\ = 1. 



,((>) 



(3) 



In Ref. [44] it was shown that Bj^ exhibits the same momentum dependence as that of the term BL , which has 
a high correlation with B^ m (about 0.936). Therefore, we further restrict our analysis only to the new non-trivial 
terms B n and B^ . In Fig. [T]we plot the shapes of the bispectrum contributions B\^ and B n . 

Let us quantify how much the shape functions of Bj^ can be fitted with the function K^ UI . We follow the procedure 
shown in Ref. [60| and calculate the shape correlator C(Bn,B' n ) defined as 



C(B u ,B' n ) 



ABn,B' n ) 



^l{B u ,B n )l{B' n ,B' n )' 



(92) 



where 



l(B n ,B' n ) = J dVMki,k2,k 3 )B' n (k u k 2 ,k 3 



h A h A h A 



(fcl 



h) 3 



(93) 



and the region of integration is defined by the condition < fci < oo, < fe/^i < 1, and 1 — kz/ki < k%/ki < 1. In 
terms of the variables kijk\ and k 3 /k\ the integral in k\ factorizes out. After performing numerical integrations we 
find that 



c(4 7) ,sr a ) 

?(7) 



-0.99989 



and 



C(4 8) , 



^k" 1 ) 



-0.99999. 



(94) 



,(«) 



Since both the contributions Bj^ and B\^ are almost completely anti-correlated with B^ , it makes sense to use 



cquil 



(7) 



3 (8) 



/ N l to measure the shapes of non-Gaussianities in the presence of the contributions B\J and B\/ as well. We also 
checked that C (B ( £\ B c £ uil ) = 0.936177, which agrees with the result in Ref. gj. 
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VI. THE NON-LINEAR PARAMETER UNDER THE SLOW- VARIATION APPROXIMATION 

In this section we derive the approximate expression of the equilateral non- linear parameter / N ^ with k\ = k,2 = ks 
on the quasi de Sitter background. In order to obtain a compact form of f^ we define the following quantities 



_ wi{Awiw 3 + 9^2) 
F 2 



^wr^x^; (95) 



12M p 4 , 



L 3X 2 P XX + 2X i P xxx + 3H0(XG 3 , X + bX 2 G 3 , xx + 2X 6 G 3 , XX x) - 2(2X 2 G 3 ^ X + X A G 3 ^ XX ) 

+6H (9X G^ xx + 16X Gi^ xxx + AX Gi tXXXX ) — SH^SXG^x + ^2X G^xx + AX Gi^ xxx ) 

+H 3 <j>(3XG 5 , x + 21X 2 G^ XX + 2AX 3 G 5 , XXX + AX A G^ XXXX ) 

-6H 2 (6X 2 G 5 ^ X + 9X 3 G 5 ^ XX + 2X A G^ XXX )] , (96) 

which are the generalizations of those introduced in Refs. |22l . 144 [45|. Notice that £ is related with Q via Q = 
4Mp 1 £/itf§. We also introduce the following quantities 

> _ XG 3 . XX _ XG^ XXX _ XG^^xxx tnr7S 

a 3X = — , Mx = — , A 5 x = — - ■ [yi) 

<^3,X <~'4,XX <-*5,XX 

We derive the leading-order terms to each /nl" coming from the coefficients Ci (i = 1, • ■ • ,8). We expand each 
/nl" m t erms °f the slow-variation parameters defined in Eq. (JTDJ), by treating c 2 , A/E, and Xi X (i = 3,4,5) as 
arbitrary parameters. Then the leading-order contributions to /^l" are given by 



/nl = -g ( ! - ^2 ) + ^2 ( e - - ^ - 4(5 G3Jf - 12<W - 32<5 G4XX + 12^6* - 10<5 G5 x - 8<5 G5 xx) , (98) 



f cquil(l) _ 10 /. 1 \ . 10 

/Nl"" (2) = ^ (^ - l) + ^ k + Vs 2s + 4fe + 2«5 G5X - 4<5 G50 ) , (99) 

/ni!" = sT( _ 2 _1 ) _ «Ty + ^7~2 (^ G3X + ^ax + 3Sg5X - 5gh - 4<5 G 50 + 8S G4XX + 25 G5XX ) 
- ^7 (3 + 2A 3X ) <5 G3 x - — (5 + 2A 4 x) S GiXX - —(4 + A 5 x)<5g5XX 

ol ol 51 

5 10 c 2 

+ ^(fc 4 + 8<5 G50 - 86 GiX - 96 G5X ) - — S. [(1 + A 3 x)<5 G3X + ^ 2 )} , (100) 

,equil(4) _ 10 £ s f-\m\ 

/NL - 27c 2' l iUi J 

/ NL Uil(5) = -j^(ts-A5 G3X -8S G4XX + 86 G5X + A6 G5XX )e s , (102) 

/ NL Uil(6) = -^~ K 1 + A 3 x )<W + 4(3 + 2A 4X )(Wx + <W + (5 + 2A 5 x)<Wx] , (103) 

ote s 

cquil(7) 65 

STL 

/ N f 1(8) = -i5§^(' s <^ f4 ' s "'^'> ) ' ( |ll5i 



/nl 1 = T7?7^ (^^ + 6<5g4X * + fcs * + 5g5 **) ' ( 104 ) 

85 



The explicit expression of the second-order term £(<5 2 ) in /S"' is given in Appendix. The above results reproduce 
those obtained for the theories with G 4 = and G5 = |45l |. 

Summing up all the terms fm" [i — 1, • • • ,8) and taking the largest contribution, we obtain 

f quii^_85A 1A 10 A 55 e s 5 y s 85 s /20 1 + A 3 x 65 

/nl 324^ c 2 y 81 E 36 c 2 12 c 2 54 c 2 ^81 e s 162c 2 e ] " G: ' ,s 

80 3 + 2A 4X 65 \ / 20 65 \ /205 + 2Asx 65 , , 

81 ^^T~ + 27^1 J fc4XX + {8U S + T62^T S ) 6g5X + {8l^7T~ + 162^1 ' Sg5XX " (1 ° 6) 
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Here we have ignored the contributions such as S G4X , 8g4<j> and Sg5<j> relative to the terms S G 4 XX /e s , S G 5 X /e s , and 

Sg5Xx/£ s - 

From Eq. (|106[1 we find that the scalar propagation speed c s mainly determines the level of non-Gaussianitics. 
Expanding the term Q in Eq. (126p in terms of slow- variation parameters, the leading-order contribution is 

Q ~ M 2 1 F[5p X (l + 2Xp X )+65 G3X {l + X 3X )-25 3 ^ + 66 G4X +245 G4XX (2 + X 4X ) 

+68 G5 x + 2S G5XX (7 + 2X 5X ) - 6S G54> ] , (107) 

where \ x (i = 3, 4, 5) are defined in Eq. ([77)) and 

Xpx = ^§^ ■ (108) 

Then the scalar propagation speed squared c 2 = M%Fe s /Q is approximately given by 

<? s ^ (Sp X + -±6 G3X - 2S G3ct> + Q8 GiX + 20S G4XX + 4£ G5 x + -±6 G5XX - 6S G5cj> )/[Sp X (l + 2\ PX ) + 6S G3X (1 + \ 3X ) 
-2S G3( k + 6S G4X + 2A5 GAXX {2 + X iX ) + 6S G5X + 2S G5XX {7 + 2X 5X ) - 6<5 G50 ] , (109) 

where we have used Eq. (|5H)) . One can estimate the ratio A/E in Eq. (|106p as follows 

— ~ [8 PX (3\p X + 2X PXX ) + 35 G3X (1 + 5X 3X + 2X 3XX ) + 6S G4XX (9 + 16A 4 x + 4:X 4XX ) + 35 G5X 

+S G5XX (27 + 24A 5 x + 4A 5XX )]/3[(5px(1 + 2A PX ) + 66 G3X (l + X 3X ) - 2<5 G30 + 6S G4X 

+245 G4XX (2 + X iX ) + 6d G5X - 6<5 G50 + 2S G5XX (7 + 2X 5X )] , (110) 



where 



v2/~< y2/°( v2/^> 

-A (JTZ.XXX . -A <^4.XXXX . -* <-t5 



, -A f^XXX . -A <^3,J>fX:X . -A ^J4.XXXX . -A (-T5.XXXX /•-i-i-n 

ApxA' — J5 , A 3 n — — , Mxx — ~ , A 5 jcx — ^ • I 111 ; 

",X t*3,X <~'4,XX <-*5,XX 

Let us consider the case in which either of the following conditions is satisfied: 

X PX > 1 , X 3X > 1 , X iX > 1 , A 5X > 1 ■ (112) 

Then one can realize c^<l and hence 1/^1"" 1^1- From Eq. (|110[) it is also possible to have A/E ^ 1 if either of 
Xpxx, X 3X x, A 4 xx, X 5X x is much larger than 1. 

VII. APPLICATIONS TO CONCRETE MODELS OF INFLATION 

In this section we apply our formula for the equilateral non-linear parameter J^l" 1° concrete models of inflation- 
such as (A) k-inflation, (B) k-inflation in the presence of the terms Gi{X) (i = 3,4, 5), (C) potential-driven Galileon 
inflation, (D) non-minimal coupling models, and (E) potential-driven inflation in the presence of the Gauss-Bonnet 
term. 

A. k-inflation 

In k-inflation in which G 3 = G 4 = G 5 = one has c 2 = 1/(1 + 2A PX ) and A/E = (1 - c 2 s )/2 + 2X PX xc 2 s /3. From 
Eq. JTP5J1 it follows that 

,equii 5 / 1 \ n ~ , „ 2n 20 a pxx , 55 e s 5 ?y s 85 s 

Since e = <5 PX = P^xX/iSM^H 2 ) [HIH, inflation occurs either around P )X -OorlsiO. 

In the former case (P,x ~ 0), as long as P. xx does not vanish at P x = 0, we have Apx = XP XX / P. x 3> 1 and 
c 2 -C 1. Hence large non-Gaussianities can be realized for the Lagrangian having a non-linear term in X. The ghost 
condensate model P = —X + X 2 /(2M 4 ) [18] belongs to this class. For the function P depending only on X there is a 
de Sitter solution at P x = 0. However this is problematic because the scalar power spectrum V-r. diverges at the de 
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Sitter solution (because c s — 0) [15J . This problem can be avoided either by involving the (^-dependence in P (such as 
the dilatonic ghost condensate model P = —X + e x ^^ Mpl X 2 / (2M 4 ) pj|) or by taking into account the terms Gi(X) 
(i = 3, 4, 5) in addition to the Lagrangian P(X) [2(|. We shall discuss the latter case in Sec. IVIIBI 

If inflation occurs in the region X sw 0, whether large non-Gaussianities can be realized or not depends on the models. 
In conventional inflation driven by the potential energy V((j>) of the field <f>, i.e. P = X — V{<j>), we have c 2 = 1 and 
Apxx = 0, which leads to the small non-linear parameter /^ ul = 55e s /36 + br) s /12. In the DBI model where the 
Lagrangian is given by P = -y/1 - 2f{(j))X/f{(j)) + 1 //(<£) - V(<j>) [3, it follows that X PX = f(4>) x /[l - 2f(4>)X] 
and e = X/[M^H 2 y/l — 2f(<j>)X]. The non-Gaussianities can be large around the region 2f(<f>)X ps 1. If the total 
energy density is dominated by the potential energy V, i.e. M 2 X H 2 ss V/3, it is possible to satisfy the condition e«l 

for X <C V (even if 2f((f>)X is close to 1). In fact this corresponds to the regime where the non-linear parameter f^ 
of the order of 10 can be achieved |17j . 

B. k-inflation with the terms Gi(X) (i — 3,4,5) 

In k-inflation models where the Lagrangian P is a function of X only, the problem of the de Sitter solution mentioned 
in Sec. IVIIAl can be circumvented by taking into account the terms Gi(X) (i = 3,4,5). For the ghost condensate 
model P = —X + X 2 /(2M 4 ) with the term G$ = [iXjM 4 (fj, > 0, M > 0) [2(J], for example, the scalar power spectrum 
is not divergent because c 2 7^ 0. In this case, in the region where the variable x — X/M 4 is close to 1, we find that 
1 - x ~ y/Sfi/MjA and /^ uil ~ 5/[6(l - a;)] ~ 4.62r~ 2 / 3 , where r ~ 16^6(1 - xf/ 2 /3. These results match with those 
found in Ref. [44] (in which detailed calculations are given). 

Let us consider the following model 



where /1 and M are positive constants. There is a de Sitter solution characterized by the condition e = Spx +6Sg4X + 
l25 GAX x = 0, at which P + 3M 2 X H 2 F - 12H 2 XG^ X = from Eq. ©. We then obtain 

R 2 = M^l~x /zM_ l-x 

36/i x ' M p 2 ! 6a; 2 (3 - 2x) ' [ ' 

where x = X/M 4 (> 0). Provided that inflation occurs in the regime \iM <C M%, the variable x is close to 1 
(with x < 1). In the following we replace x for 1 apart from the terms including 1 — x. From Eq. (|115p we have 
fiM/M^ ~ (1 - x)/6 and H 2 ~ M 4 /(6M 2 {). Since Q ~ 12M 2 { > 0, the no-ghost condition is satisfied. The scalar 
propagation speed squared is given by c 2 ~ 2(1 — x)/9, so that the Laplacian instability is absent for x < 1. Since 
e s ~ 8(1 — x)/3, the scalar power spectrum and the tensor-to-scalar ratio are given by 

p *>/*_( M_X 1 r ^i^(l-,)3/ 2 (116) 

' n ~ 2567T 2 [m p J (1-x)3/2' r - 9 ( } ■ [ ' 

The term A/E in Eq. (|110p is estimated to be A/E ~ 1/2. Since Sqaxx — (1 — x)/3, we find that the equilateral 
non-linear parameter is given by 

25 1 
144 ~x 



f^~J2 7T ±-~ 1.28r- 2 / 3 . (117) 



For smaller r, /nl" gets larger. If r = 0.01, then /^l" — 9-4- This level of non-Gaussianities can be detectable in 
future observations. 



P=- X +—,, G 5 = M7 — . (118) 



We also consider the following model 

X 2 X 2 

2M 4 ' 5_/i Mi 

A similar calculation shows that there is a de Sitter solution characterized by jjtM /M% ~ 27(1 — x) 2 /25 and 
H 2 ~ Af 4 /(6M 2 1 ) for a; = X/M 4 close to 1. Since e s ~ 18(1 - x)/5, c 2 ~ 3(1 - x)/10, £ G 5X = fcsxx ^ 6(1 - x)/5, 
A/E ~ 1/2, and Q ~ 12Af 2 j > in this case, it follows that 

f cc i uil ^ * 17r~ 2/3 (119) 

/nl " 1458 l-x ' ( ' 

which is about one order of magnitude smaller than (|117[) . 
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C. Potential-driven Galileon inflation 

Let us proceed to the potential-driven inflation [P = X — V((f>)] in the presence of the Galileon term Gs(X) ex X n 
(n > 0) [33| . The inflationary dynamics for the case n = 1 were studied in Refs. |33l . |38j. Since \px — and 
A3 x = n — 1 we have 

c a = <W+4te (120) 

opx + bnd G3 x 

In the regime Sq3X 3> 5px it follows that c 2 ~ 2/(3n) and hence c 2 <S 1 for n ^> 1. In this case one has e s ~ 4<5g3X 
and A/S ~ n/3, so that the non-linear parameter f|106[) is given by 

/Nr--S|- (121) 



In the presence of the term G±(X) oc X n the slow- variation parameter Sqaxx is related with <5g4X, via <5g4XX = 
(n — \)5gax- If n 3> 1, one has c^ ~ 5/(6n), A/S ~ n/3, e s ~ 20<5g4xx in the regime 1 ^> Sgaxx 3> <5px, and hence 
/n1 u11 ^ -137n/1215. Similarly in the case where G 5 (X) ex AT" with n > 1, it follows that /^ uil ~ -155n/1944 for 

1 > £ G 5XX > Sp X - 

D. Nonminimal coupling models 

The nonminimal coupling of the scalar field <f> with the Ricci scalar R corresponds to the choice G4 = F{4>), where 
F{4>) is an arbitrary function in terms of <f>. In the absence of the terms G3 and G5, using <5g4X = fc«x = 0, the 
scalar propagation speed squared is c 2 = 1/(1 + 2\px)- Hence, as in the case of k-inflation, we require Xpx 3> 1 to 
realize large non-Gaussianities. For the theories in which P does not have non-linear terms in X we have c~ = 1 and 

f^t l = 0(e s , Vs ). (122) 

For example, Brans-Dicke theories described by the action P — usbdM p \X/(/) — V{(j)) [48] belong to this class. Hence 
the non-Gaussianities in those theories are small (including f(R) gravity where the Brans-Dicke parameter ojbd is 

0). 

The models with field derivative couplings recently studied in Refs. [49l, l50fl correspond to G5 = F (</>), in which case 
$G5X — and <5g5XX = 0. In the absence of the terms G3 and G4 we have 

c 2 = §px - 6^ G5 ^ 123 ^ 

SpxO- + 2A PX ) - 66054, 

If the Lagrangian does not include non-linear terms in X (like the models discussed in Ref. [49]), then c 2 = 1 and 
/nl" = C( e s7 7 7s)- This is consistent with the results found in Ref. [50]. 

E. Potential-driven inflation with a Gauss-Bonnet term 

The action ([T]) even covers the Gauss-Bonnet coupling of the form —£(4>)G, where Q = R 2 — 4R a pR al3 + R a /j 1 sR a ^ lS 
is the Gauss-Bonnet term. If one chooses the following combination of P, G3, G4, and G5 

P = -8£ (4) (<W^ 2 (3 - lnX) , G 3 = -A^\(f)X{l - 31nX) , 

G 4 = -4£,^(<j>)X{2-\nX), G 5 = 4^ (1) (0) hxX , (124) 

where £^(0) = d n £((f>)/d(j> n , one can show that the field equations following from this Lagrangian are equivalent to 
those derived by the Lagrangian —£(4>)G [43]. 

Let us consider the case of potential-driven inflation in which the Lagrangian X — V(4>) is added to P in Eq. (|124|l . 
Using the choice of the functions in Eq. (|124[) , we find that the scalar propagation speed squared is given by 

MS 2 
c 2 s = l-—^(65 i + S x ), (125) 

Ox 
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where Sx = X/(M 2 l H 2 ) and <% = H^/M 2 ^ This expression agrees with the one found in Ref. j38|. Since c 2 — 1 is a 
second-order quantity, the scalar propagation speed is very close to 1. Furthermore we can show that, in this case, 
A/S = 0, and that the terms in Eq. (|106[) coming from the functions S G3 x, S G4 xx, Sq5X, and 5q5XX are of the order 
of e, namely proportional to 5^ (4e + 2r]£ — rjx)/Sx- Using e s ~ 5x, the leading contribution to /JS" 1 is given by 

/N q L Uil = §S X + ^vx + ^(4e + 2 Vi - r, x ) , (126) 

where r/£ — S^/(HS^) and r\x — 5x/{H5x)- Hence the non-linear parameter in this model is small. 

VIII. CONCLUSIONS 

In this paper we have evaluated the primordial non-Gaussianities generated during inflation for the models described 
by the action ((T]). Our analysis is general enough in that it covers a wide variety of single-field models having second- 
order equations of motion. 

The procedure to obtain the three-point correlation function of curvature perturbations is similar to that given in 
Ref. [45| . There are 8 shape functions which contribute to the scalar non-Gaussianities. Among them, five functions are 
already known to give rise to the equilateral type of non-Gaussianities. We studied the shapes of the remaining three 
functions and found that they can be well approximated by the same type as well. Hence the dominant contribution 
to the three-point correlation function comes from the equilateral configuration. 

We derived the equilateral non-linear parameter /nl" un der the approximation that the slow-variation parameters 
defined in Eq. (JTUJ) are much smaller than 1. The formula (|106[) is valid for any quasi de Sitter background and thus it 
is convenient to apply it to concrete single-field models of inflation. In fact we used this formula to a number of models 
such as (A) k-inflation, (B) k-inflation with the terms Gi(X) (i = 3,4,5), (C) potential-driven Galileon inflation, (D) 
nonminimal coupling models without a non- linear term in X, (E) potential-driven inflation with a Gauss-Bonnet term. 
In the models (D) and (E) we showed that c 2 is close to 1 and hence the non-Gaussianities are small. However, for the 
models (A), (B), (C), it is possible to realize I/^l" I ^" 1 depending on the choice of the functions P,Gi (i — 3,4,5). 

The potential detectability of non-Gaussianities in future observations will allow us to distinguish between different 
inflationary models. In particular we expect that the precise measurement of /nl" as well as n-ji and r will provide 
significant implications for the viability of single-field models in which c 2 is much smaller than 1 . 
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Note added- While we were completing this work, we became aware of the paper by Gao and Steer [6J| who calculated 
primordial non-Gaussianities in the same model as ours. Their formula for the bispectrum (98) is consistent with our 
formula (155)) . 

Appendix A: The second-order term in f^ 1 

The second-order term £(<5 2 ) in Eq. poop is given by 

£,(S 2 ) = [(6 + 4X 3X )5 G 4x + 8(3 + X 3X + A 4X )<Wx + 3(2 + X 3 x)S G 5X + (9 + 2A 3 x + "2X 5X )5 G5 xx - fc^(l + hx) 
-2(5 G50 (3 + 2X 3X )}S G 3X + 8S 2 G4X + [16(5 + 2X 4X )S G 4XX + 18S G5X + 8(4 + X 5X )S G5X x - 2S G4<f> 
-165 G ^}5 G4X + 64(2 + X 4X )S 2 G4X x + [24(3 + X 4X )<W + 8(11 + 2X 4X + 2X 5X )S G5X x ~ 8<W2 + X 4X ) 
-16<5 G50 (5 + 2X 4 x)]8 G4 xx + 9S 2 G5X + [3(9 + 2X 5X )8 GbX x - 35 GH - 18S G5 ^S G5X + 2(7 + 2X 5X )S 2 G5XX 
-[(7 + 2X 5X )S G44> + 8(4 + A 5 x)<5 G 5<#G5XX + 8S 2 G5<f> + 2S G4 ^S G5<P . (Al) 
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